March 8, 2013 0:36 World Scientific Review Volume - 9in x 6in bcs50 



Chapter 1 
Microscopic origin of pairing 



Eduard E. Saperstein and Marcello Baldo 

Kurchatov Institute, Moscow; saper@mbslab.kiae.ru 
INFN Catania; baldo@ct.infn.it 

A brief review of recent progress in the ab intio theory of nuclear pairing 
is given. Nowdays several successful solutions of the ab intio BCS theory 
gap equation were published which show that it is a promising first step 
in the problem. However, the role of many-body correlations that go 
beyond the BCS scheme remains uncertain and requires further investi- 
gations. As an alternative, the semi-microscopic model is discussed in 
which the effective pairing interaction calculated from the first princi- 
ples is supplemented with a small phenomenological addendum contain- 
ing one phenomenological parameter universal for all medium and heavy 
atomic nuclei. 



1. Introduction 

Recently, the fifty years anniversary of Cooper pairing in nuclei 1 took place. 
However, only in the last few years some progress has been made in the 
microscopic theory of nuclear pairing, first, by the Milan group, 2-4 a little 
later by Duguet et al. 5 ' 6 and finally by the Moscow-Catania group. 7-11 
In the first seminal paper of the Milan series, the BCS gap equation for 
neutrons with the Argonne V14 potential was solved for the nucleus 120 Sn. 
The Saxon- Woods Shell-Model basis with the bare neutron mass m* = m 
was used, and the discretization method in a spherical box was applied to 
simulate the continuum states restricted by the limiting energy -E max = 600 
MeV. Rather optimistic result was obtained for the gap value, Abcs = 2.2 
MeV. Although it is bigger of the experimental one, A exp ~ 1.3 MeV, the 
difference is not so dramatic and left the hope to achieve a good agreement 
by developing corrections to the scheme. In Refs. 3,4 the basis was enlarged 
to -Emax = 800 MeV, and, what is more important, the effective mass 
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to* 7^ m was introduced into the gap equation. The new basis was calcu- 
lated within the Skyrme-Hartree-Fock (SHF) method with the Sly4 force, 12 
that makes the effective mass m*(r) coordinate dependent and essentially 
different from the bare one to. E.g., in nuclear matter the Sly4 effective 
mass is equal to to* = 0.7m. As it is well known, in the weak coupling limit 
of the BCS theory, the gap is exponentially dependent, i.e. A oc exp(l/g), 
on the inverse dimensionless pairing strength g — m*V e ffkF /w 2 , where V e ff 
is the effective pairing interaction. Therefore, a strong suppression of the 
gap takes place in the case of m* < m. The value of Abcs = 0.7 MeV 
was obtained in Ref. 3 and Abcs = 1-04 MeV, in Rcf. 4 In both cases, the 
too small value of the gap was explained by invoking various many-body 
corrections to the BCS approximation. The main correction is due to the 
exchange of low- lying surface vibrations ("phonons"), contributing to the 
gap around 0.7 MeV, 3 so that the sum turns out to be A = 1.4 MeV very 
close to the experimental value. In Ref., 4 the contribution of the induced 
interaction caused by exchange of the high-lying in- volume excitations was 
added either, the sum again is equal to A ~ 1.4 MeV. Thus, the calcu- 
lations of Refs. 3 ' 4 showed that the effects of to* ^ to and of many-body 
corrections to the BCS theory are necessary to explain the difference of 
(Abcs — A oxp ). In addition, they are of different sign and partially com- 
pensate each other. Unfortunately, both effects contain large uncertainties. 
This point is discussed in Section 3. 

In 2008 Duguet and Losinsky 5 made a fresh insight to the problem 
by solving the ab initio BCS gap equation for a lot of nuclei on the same 
footing. It should be noticed that the main difficulty of the direct method to 
solve the nuclear pairing problem comes from the rather slow convergence 
of the sums over intermediate states A in the gap equation, because of 
the short-range of the free TV TV-force. To avoid the slow convergence, the 
authors of Refs. 5 ' 6 used the "low-k" force Vi ow -k 13 ' 14 which is in fact very 
soft. It is defined in such a way that it describes correctly the AA-scattcring 
phase shifts at momenta fc<A, where A is a parameter corresponding to the 
limiting energy ~ 300 MeV. The force Vi ow -k vanishes for fc>A, so that in 
the gap equation one can restrict the energy range to -E max ~300 MeV. In 
addition, a separable version of this force was constructed that made it 
possible to calculate neutron and proton pairing gaps for a lot of nuclei. 
Usually the low-k force is found starting from some realistic AA-potential 
V with the help of the Renormalization Group method, and the result does 
not practically depend on the particular choice of V. 13 In addition, in Ref. 5 
Viow-k was found starting from the Argonne potential vi 8 , that is different 
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only a little bit from Argonne V14, used in Ref. 4 Finally, in Ref. 5 the same 
SLy4 self-consistent basis was used as in Ref. 4 Thus, the inputs of the two 
calculations look very similar, but the results turned out to be strongly 
different. In fact, in Ref. 5 the value Abcs — 1-6 MeV was obtained for the 
same nucleus 120 Sn which is already bigger than the experimental one by 
~ 0.3 MeV. In Refs. 7,8 the reasons of these contradictions were analyzed. 
It turned out that, in fact, these two calculations differ in the way they 
take into account the effective mass. It implies that the gap A depends not 
only on the value of the effective mass at the Fermi surface, as it follows 
from the above exponential formula for the gap, but also on the behavior of 
the function m*(k) in a wide momentum range. However, this quantity is 
not known sufficiently well, 8 which makes rather uncertain the predictions 
of such calculations. 

To avoid such uncertainties, a semi-microscopic model for nuclear pair- 
ing was suggested by the Moscow-Catania group. 9-11 It starts from the 
ab initio BCS gap equation with the Argonne force vi§ treated with the 
two-step method. The complete Hilbert space S of the problem is split into 
the model subspace So of low-energy states and the complementary one S'. 
The gap equation is solved in the model space with the effective interac- 
tion V e ff which is found by projecting out the complementary subspace. A 
new version of the local approximation, the so-called Local Potential Ap- 
proximation (LPA), 15 is used in the subspace S'. This ab-initio term of 
Vefi is supplemented by a small addendum proportional to the phenomcno- 
logical parameter 7 that should hopefully embody all corrections to the 
simplest BCS scheme with m* = m. Smallness of the correction term is 
demonstrated in Fig. 1 where a localized "Fermi average" form of V ff is 
displayed without (7 = 0) and with (7 = 0.06) the phenomenological cor- 
rection. Non-negligible effect of so small change of V e ff to the gap value 
is owing to the above mentioned exponential enhancement effect. Explicit 
definition of the functions displayed in Fig. I is given in Section 4 where 
some results of the semi-microscopic model are presented. 

Fig. 1 demonstrates also the surface nature of nuclear pairing, the 
effective pairing interaction at the surface is ten times stronger than inside. 
This explains the relative success of the BCS approximation. Indeed, at 
the surface the main corrections to the Brueckner-like theory which leads 
to the BCS scheme are small. 
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Fig. 1. The Fermi average effective pairing interaction V ( ^j(i?) for 120 Sn and 2ut) Pb 
nuclei 

2. The microscopic BCS equation. LPA approximation 

The general many-body form of the equation for the pairing gap is as fol- 
lows, 16 

A T =U T G T G S T A T , (1) 

where r = (n,p) is the isotopic index, IA T is the AW-interaction block 
irreducible in the two-particle r-channel, and G T {G S T ) is the one-particle 
Green function without (with) pairing. A symbolic multiplication denotes 
the integration over energy and intermediate coordinates and summation 
over spin variables as well. The BCS approximation in Eq. (JTJ means, first, 
the change of the block U of irreducible interaction diagrams with the free 
AW-potential V in Eq. (p}, and, second, the use of simple quasi-particle 
Green functions G and G s , i.e. those without phonon corrections and so 
on. In this case, Eq. ([T]) turns greatly simplified and can be reduced to the 
form usual in the Bogolyubov method, 

A r = -V T x T , (2) 

where 



x T = J^G T G°A T 



(3) 
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is the anomalous density matrix which can be expressed explicitly in terms 
of the Bogolyubov functions u and v, 



Summation in Eq. Q scans the complete set of Bogolyubov functions with 
eigen-energies Ei > 0. 

To overcome the slow convergence problem, a two-step renormalization 
method of solving the gap equation in nuclei was used in Refs. 9 ~ n The com- 
plete Hilbert space of the pairing problem S is split in the model subspace 
Sq, including the single-particle states with energies less than a separation 
energy Eq, and the complementary one, S' . The gap equation is solved in 
the model space: 



with the effective pairing interaction V^ s instead of the block V T in the BCS 
version of the original gap equation (fTJ). It obeys the Bethe-Goldstone type 
equation in the subsidiary space, 



In this equation, the pairing effects can be neglected provided the model 
space is sufficiently large, Eq ;§> A. That is why we replaced the Green func- 
tion G S T for the superfluid system with its counterpart G T for the normal sys- 
tem. To solve Eq. (j6|) in non-homogeneous systems a new form of the local 
approximation, the Local Potential Approximation (LPA), was developed 
by the Moscow-Catania group. Originally, it was found for semi-infinite 
nuclear matter, 17 then for the slab of nuclear matter (see review article 15 ) 
and, finally, for finite nuclei. 7:8 It turned out that, with a very high accu- 
racy, at each value of the average cm. coordinate R = (rj + r 2 + r 3 + r^/A, 
one can use in Eq. ([6]) the formulae of the infinite system embedded into the 
constant potential well U = U (R) . This significantly simplifies the equa- 
tion for V c ff , in comparison with the initial equation for A. As a result, the 
subspace S' can be chosen as large as necessary to achieve the convergence. 
Accuracy of LPA depends on the separation energy Eq. For finite nuclei, 
the value of Eq=4Q MeV guarantees the accuracy better than 0.01 MeV for 
the gap A. 

Let us notice that the use of the low-k force Vi ow _k could be also inter- 
preted in terms of the two-step renormalization scheme of solving the BCS 
version (U— >V) of the gap equation (fTJ), with i?o— 300 MeV and with free 
nucleon Green functions G in Eq. © (i.e. U(R) = 0). Then, one obtains 




(4) 



A r = V^G T G 6 T A T \ Sl 



(5) 



Vl s = V T + V T G r G T Vl s \ s ,. 



(6) 
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Voff^Vjow-k (see Ref. 14 where the usual renormalization scheme is used to 
find V|ow-k instead of the Renormalization Group equation). 

3. Corrections to the BCS scheme 

As it was mentioned in the Introduction, there are mainly three types of 
corrections to the plain ab initio BCS gap equation with bare nucleon mass 
m. The first one is the effect of the effective mass m* ^ m considered in 
Refs. 3 ' 4 and Refs. 5,6 as well. The second one is the contributions of low- 
lying surface phonons 3,4 and the third one, the induced pairing interaction 
due to high-lying in-volume excitations. 4 



Fig. 2. Momentum dependence of the effective mass m*(k) at different Fermi momen- 
tum values kp. 

Let us begin from the effective mass. In the analysis 7,8 of the difference 
between the BCS gap values of Ref. 4 and Refs., 5,6 it was found that the 
gap A depends on the behavior of the function m*(k) in a wide momentum 
range. However, this quantity is not known sufficiently well even in nuclear 
matter. In Fig. 2 the effective mass m*{k) of symmetric nuclear matter 
at different Fermi momentum values kv are displayed. They were found 8 
with the Brueckner-Hartree-Fock method which is also, of course, an ap- 
proximation. These functions behave rather non-regular in vicinity of kp 
and tend to the bare mass value very slowly. This makes rather doubtful 
the way to account for the m* effect in Refs. 3,4 and in Refs. 5 ' 6 as well. In 
the first case, the effective mass was taken fc-independent (equal to Wg Ly4 ) 
till the cut-off fc max ~6 fm _1 . In the second case, the ansatz was used of 
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m*=77ig Ly4 for fc<A~3 fm -1 and m*=m for fc>A. 

Similar problems appear if one tries to use an explicit form of the Z- 
factor, which is an additional ingredient of the gap equation, not only in the 
combination yielding the complete effective mass m*(k). 8 In this Section, 
speaking for brevity about the effective mass we mean in fact both the 
fc-mass and .E-mass. 

Corrections to the BCS scheme owing to phonon contributions are of 
primary importance. To our knowledge, the most advanced calculation 
of these corrections was carried out in Ref. 3 It includes as the pho non- 
exchange term of the block U in Eq. ([1]), the so-called induced interaction, 
and corrections to the Green functions G, G s as well. However, it also 
possesses some deficiency connected with omitting the so-called tadpole 
diagrams. Up to now, they were consistently taken into account within 
the self-consistent Finite Fermi Systems (FFS) theory 18 only for magic 
nuclei where pairing is absent. It turned out that their contribution is 
usually important and often is of the opposite sign to the usual diagrams 
diminishing the total value of the effect under consideration. This formalism 
was generalized for superfluid nuclei in Ref., 19 but numerical applications 
are still absent. 

Let us finally discuss corrections to the BCS version of Eq. (JlJ due to 
the induced interaction from high-lying particle-hole in-volume excitations. 
The attempt in Ref. 4 to determine the latter from the SLy4 force together 
with the nuclear mean field looks questionable. Indeed, the SLy4 parame- 
ters were fitted to the nuclear mass table data mainly related to the scalar 
Landau-Migdal (LM) amplitudes /, /'. As to the spin amplitudes g,g', 
they remain practically undetermined in the SHF method. But the contri- 
bution of the spin channel to the induced interaction is not smaller than 
that of the scalar one. 4 The LM parameters g, g 1 are well known from the 
calculations of nuclear magnetic moments within the FFS theory 20 but, as 
for the Skyrme parameters, only at the Fermi surface. However, the states 
distant from the Fermi surface are important to calculate the induced inter- 
action. The induced interaction for such states has been determined only 
in nuclear matter within the microscopic Brucckner theory. 21 

4. The semi-microscopic model for nuclear pairing 

To avoid uncertainties of explicit consideration of corrections to the BCS 
scheme discussed above, the semi-microscopic model was suggested in 
Refs. 9 ~ n In this model, a small phenomenological addendum to the ef- 
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fective pairing interaction is introduceed which embodies approximately all 
these corrections. The simplest ansatz for it is as follows: 

VJ ff (n, r 2 , r 3 , r 4 ) = 7 T B e f (n, r 2 , r 3 , r 4 ) + l T c/-^- f[ c5(n - r<). (7) 

Here p(r) is the density of nucleons of the kind under consideration, and 7 1 " 
are dimensionless phenomenological parameters. To avoid any influence of 
the shell fluctuations in the value of p(0), the average central density p(0) 
is used in the denominator of the additional term. It is averaged over the 
interval of r<2 fm. The first, ab initio, term in the r.h.s. of Eq. ([7]) is the 
solution of Eq. © in the framework of the LPA method described above, 
with m*—m in the subspace 5". 
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Fig. 3. Neutron gap in Sn isotopes 



We will see below that a rather small value of the phenomenological 
parameter j n — 7 P ~ 0.06 is sufficient to produce the necessary effect of 
suppressing theoretical gaps predicted by the ab initio calculation. The 
smallness of the phenomenological addendum to the effective interaction 
itself is demonstrated in Fig. 1 where the localized "Fermi average" effective 
interaction is drawn for 7 = and 7 = 0.06 values for two heavy nuclei. In 
the mixed coordinate-momentum representation, this quantity is defined as 
follows: V c ff(ki,k 2 ,ri,r 2 ) -> V^ S (R = n)H i=2 S(ri - i-*), where 



V e F ff (i?)= d 3 tV ee (k 1 = k 2 = k F (R),R-t/2,R + t/2), (8) 
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with k F (R) = ^2m{n - U(R)), provided fi - U{R) > 0, and k F (R) = 
otherwise. Here \i and U(R) are the chemical potential and the potential 
well, respectively, of the kind of nucleons under consideration. A similar 
quantity was considered before in the nuclear slab to visualize the effective 
interaction properties. 15 
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Fig. 4. Neutron gap in Pb isotopes. 



In Ref., 10 the above equations were solved in the self-consistent A-basis 
of the Energy Density Functional (EDF) of Fayans et al. 22 ' 23 Two sets of 
the functional were used, the original one DF3 23 and its modification DF3- 
a. 24 In the latter, the spin-orbit and effective tensor terms of the initial 
functional were modified. The results for the pairing gap in three chains of 
semi-magic nuclei are displayed in Figs. 2-4. 

In accordance with the recipe of Ref., 4 we represent the theoretical gap 
with the "Fermi average" combination 



A F =^(2j+l)A AA /^(2. 7 +l), 



(9) 



where the summation is carried out over the states A in the interval of 
\e\— ^i|<3 MeV. The "experimental" gap is determined by the symmetric 
5-term odd-even mass difference. As it is argued in Ref., 10 the relevance 
of the mass difference to the gap has an accuracy of ~ (0.1 -j- 0.2) MeV. 
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Fig. 5. Proton gap for N = 82 isotones 



Therefore, it is reasonable to try to achieve the agreement of the gap within 
such limits. It should be noted also that the theoretical accuracy of the 
approach based on the "developed pairing" approximation 16 with particle 
number conservation on average only is also about 0.1 MeV. 25 

Let us begin from neutron pairing and consider first the tin isotopes, Fig. 
2. We see that the BCS gap (7=0) is approximately 30% greater than the 
experimental one. Switching on of the phenomenological addendum with 
7=0.06 makes theoretical gap values closer to experiment. However, predic- 
tions of two versions of the functional used are significantly different, being 
much better for the DF3-a functional. In particular, the A-dependcnce of 
the experimental gap is reproduced with a pronounced minimum in the 
center of the chain. As the analysis in Ref. 10 has shown, this strong differ- 
ence between results for two functionals is ought to the strong influence to 
the gap of the high j intruder state l/in/2- Its position depends essentially 
on the spin-orbit parameters and is noticeably different for DF3 and DF3-a 
functionals. It explains the effect under discussion. 

In the lead chain, see Fig. 3, the overall pattern is quite similar. Again 
the BCS gap is approximately 30% bigger of the experimental one and again 
inclusion of the phenomenological term with 7=0.06 gives a qualitative 
agreement. Now, the difference between two functionals is much less. In 
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this case, the agreement is quite perfect for the DF3 functional and a little 
worse for the DF3-a one, but also within limits for the accuracy discussed 
above. 

Let us go to proton pairing, N=82 chain, see Fig. 4. In this case, the 
Coulomb interaction should be included into the pairing effective interac- 
tion, 



As it is argued in Ref., 10 the bare Coulomb potential could be with high 
accuracy used in this equation. The strong Coulomb effect in the gap is 
demonstrated in Fig. 4. It is also explained with the exponential depen- 
dence of the gap on the pairing interaction. It should be mentioned that 
Duguet and co-authors 6 were the first who inserted the Coulomb interac- 
tion into the pairing force for protons. Only after inclusion of the Coulomb 
interaction into V e ff, we can use the same value of 7 = 0.06 for protons 
and neutrons. As for tin isotopes, the difference between DF3 and DF3-a 
results is rather strong, now in favor of the DF3 functional. This effect is 
again due to different positions of the l/in/2 level, but now for protons. 
Overall agreement with experiment is for protons worse, maybe, because of 
closeness of some nuclei to the region of the phase transition to deformed 
state. 

5. Conclusions 

We reviewed briefly the recent progress in the microscopic theory of pair- 
ing in nuclei involving contributions of the Milan group 2-4 of Duguet with 
coauthors 5,6 and of Moscow-Catania group. 9-11 It became clear that the 
plain ab initio BCS gap equation with bare mass m* = m 2,10 is a good 
starting point for such theory. As the analysis in Rcfs. 3 ' 4 showed, the ef- 
fect of m* 7^ to, and that of many-body corrections to the BCS theory are 
of different sign and partially compensate each other. Qualitatively, they 
can explain the difference between Abcs and A exp but both effects contain 
large uncertainties and hardly can be took into account definitely at the 
modern level of nuclear theory. 

As an alternative, the semi-microscopic model was suggested in 
Rcfs. 7 ' 9 ' 10 starting from the plain ab initio BCS theory with the use of 
the self-consistent EDF basis characterized by the bare nucleon mass. The 
gap equation is recast in the model space So, replacing the bare interaction 
with the effective pairing interaction V e s determined in the complementary 




V " ff + V C . 



(10) 
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subspace S'. The Argonne v 18 potential was adopted to find V e ff along 
with the LPA method. 15 A small phcnomenological term is added to this 
effective interaction that contains one parameter, common to neutrons and 
protons, which should embody approximately the effective mass and other 
corrections to the pure BCS theory. Calculations were carried out with two 
versions of the EDF, the initial DF3 functional 23 and its version DF3-a 24 
with modified spin-orbit and effective tensor terms. They involve semi- 
magic lead and tin isotopic chains and the N = 82 isotonic chain as well. 
The Coulomb interaction is explicitly included in the proton gap equation. 
It was found that the model reproduces reasonably well the experimental 
values of the neutron and proton gaps for both functionals. However, the re- 
sults depend essentially on the single-particle spectrum of the self-consistent 
basis used, especially on the position of high j-levels. Thus, for tin isotopes 
agreement is much better for the DF3-a functional which reproduces better 
the position of the "intruder" l/in/2 neutron level. On the contrary, for 
N = 82 isotones agreement for this functional is worse. In this case, the 
DF3 functional describes better the position of the same intruder state but 



for protons. The overall disagreement is y (A t h — A cxp ) 2 ~0.13 McV for 
the DF3 functional and ~0.14 MeV, for the DF3-a one. 

The model under discussion exhibits a week point by including to the 
"universal" phenomenological addendum all corrections to the BCS scheme. 
The effective mass correction and the one due to the induced interaction 
from high-lying excitations are mainly in-volume. Therefore, indeed, they 
should be universal for medium and heavy nuclei. On the contrary, the 
phonon correction is surface and may vary from one nucleus to another as 
low-lying phonon characteristics do. A more consistent scheme should, ev- 
idently, include the explicit consideration of the low-lying phonons, as e.g. 
in, 3 but taking into account the tadpole diagrams. 19 In this case, the phe- 
nomenological constant 7 should, of course, change. Such extended model 
is much more complicated but should be more accurate in reproducing ex- 
perimental data. 
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